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ABSTRACT 

We present optical (UBVJ), ultraviolet (FUV, NUV), and high resolution 
atomic hydrogen (H I) observations of the nearby blue compact dwarf (BCD), 
VII Zw 403. We find that VII Zw 403 has a relatively high H I mass-to-light 
ratio for a BCD. The rotation velocity is nominally 10 — 15 km s^^, but rises to 
~ 20 km s~^ after correction for the ~ 8 — 10 km s~^ random motions present 
in the gas. The velocity field is complex; including a variation in the position 
angle of the major axis going from the NE to the SW parts of the galaxy. Our 
high resolution H I maps reveal structure in the central gas, including a large, 
low-density H I depression or hole between the southern and northern halves of 
the galaxy, coincident with an unresolved x-ray source. 

Although interactions have been proposed as the triggering mechanism for the 
vigorous star formation occurring in BCDs, VII Zw 403 does not seem to have 
been tidally triggered by an external interaction, as we have found no nearby 
possible perturbers. It also doesn't appear to fall in the set of galaxies that 
exhibit a strong central mass density concentration, as its optical scale length is 
large in comparison to similar systems. However, there are some features that 
are compatible with an accretion event: optical/H I axis misalignment, a change 
in position angle of the kinematic axis, and a complex velocity field. 

Subject headings: galaxies: irregular — galaxies: star formation — galaxies: 
individual: VII Zw 403 — galaxies: ISM 
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Introduction 



VII Zw 403 is cla ssified as a bl ue compact dwarf (BCD) galaxy. Compact galaxies 
were first identified by IZwickyI (jl964f ) as high surface brightness objects that were just dis- 
tinguishable from stars on the Palomar Sky Survey. This origin al classification included 
galaxies with a wide range of colors and total luminosities. Later, iThuan &: Martini (jl98l[ ) 
presented a study of 115 blue compact dwarf galaxies, which they defined as low luminosity 
{Mb < —18) galaxies with a compact optical size (< 1 kpc in diameter), that have optical 
spectra with "strong sharp narrow emission lines superposed on a blue continuum." The 
last two properties were chosen to ensure that the objects in their study had dense regions 
with active star formation. 

It was shown early on that B CD galaxies couM not h ave maintained such vigorous star 
formation for long periods of time ( ISearle &: Sargentlll972l ). so such activity must be episodic. 
The cause of this activity is unclear. Unlike spiral galaxies, triggering mechanisms such as 
spiral density waves and shear are not present in dwarf galaxies, which tend to be slow, 
mostly solid-body rotators. Gravitational instabilities in the gas, which work well to predict 
the location of star formation in large galaxies, se em to work differe ntly in dwarfs, where the 
average gas density is typically below the critical iKennicuttI (119891 ) value by a factor of 2 or 
3, and yet, stars still form. 



Hunter fc Elmegreenl (J2004j ) showed that the nature of the H II regions in BCDs are not 
that different from those in Ims; they are just crowded more closely together towards the 
center in the BCD galaxies. The amount of star formation caused by pressurized triggering 
versus spontaneous instabilities is likely similar as well. So what causes the centralized star 
formation bursts that characterize BCDs? Two possibilities explored in this paper incl ude 
an inherent central mass density concentration (jSalzer et al.l Il999l : Ivan Zee et al.l 120011) or 



gravitational interac t ions with nearby objects, i ncluding gas infall/ accretion (JTaylorl 119971 : 
Pustilnik et aPboOll : iHunter fc Elmegreenl 12004 : iBrosch et aDl2004l ). 



Although gravitational interactions are often invoked to explain bursts of star forma- 
tion in galaxies, previous studies have shown that not all the star formation activity in 
BCDs can be explained by external gravitational triggers. A plot of current normalized 
star f ormation rate vs. distance to the nearest projected (catalogued) galaxy (see Figure 
14 in iHunter fc Elmegreenl |2004| ) shows no obvious correlation between presence/distance 
to companion and star formation activity for Ims, BCDs, or Sms in a sample of 140 galax- 



ies. This is consistent with other searches for bright c ompanions ( ICampos-Aguilar fc Moles 



1991 



Eq3) 



Telles fc Terlevichlll995t iTelles fc Maddoxll2000l ) . and with work done by 



Brosch et al. 



This does not rule out possible triggering by much fainter compani ons, though 



and there is some evidence for an excess of small/faint companions to BCDs (JTaylorl Il997l : 
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Pustilnik et al.l 120011). and there are several examples of BCD s with nearby small galaxies 
or gas clouds (e.g. IPutman et al.l llQQSt IPustilnik et al.ll2003[ ). Such faint companions are 
dif ficult to detect and may ha ve been missed in earlier studies of BCDs. As pointed out 
by iHunter fc Elmegreenl (120041 ). dwarf galaxies are common so "the most likely object to 
perturb any galaxy is a dwarf," but the effect of a small/faint companion on a nearby dwarf 
galaxy will be more dramatic than on a much larger system. They suggest that a possible 
mechanism for driving mass inflow in dwarfs would be perturbation by a small companion, 
resulting in angular momentum redistribution, and that such interactions could have hap- 
pened long ago. Once a low angular momentum system is created, it can continue to host 



episodic star formation for long time periods (e.g. Ivan Zee et al.l 1200 ll ) 



The re is some evidence for central mass density concentration in BCDs in both t he neu- 



tral gas fiTavlor et al.lll994j . Il995 : 



van Zee et al. 



19981: 



the op tical underlying host galaxy (jPapaderos et al.lll996a| ]bl: ISalzer et al.lll999l ) . IVan Zee et al 



Simpson &: GottesmanI l2000f) and 



m 



(J200ll ) find that the 6 BCDS in their sample have steep rotation curves relative to low surface 
brightness dwarfs (LSBDs) with similar luminosities, with their mass distributions (stars, 
gas, and dark matter) more centrally concentrated, and that they are lower angular momen- 
tum systems than similar LSBDs. They suggest that it is just such low angular momentum 
systems that are able to collapse into centrally concentrated objects with high central densi- 
ties. For a solid body rotation curve, the star formation threshold de nsity is propor tional to 
the slope of the rotation curve in a Toomre disk instability analysis ( lToomrelll964l ). Hence, 
a steep rotation curve raises the critical threshold density for star formation. This delays 
the onset of star formation until the gas densities are quite high; and once star formation 
does begin, the system is more susceptible to bursting behavior than a more diffuse (higher 
angular momentum) system. 

Another possible cause for the star formation activity that characterizes BCDs is the 
accretion/merger of a small gas cloud. Such events have been invoked to explain mismatches 
between the optical and (gas) kinen aatic axes in galaxies, t wists in the velocity field, ex- 
tended/unusual gas morphology (e.g. iHunter fc Wilcotsll2002l ). counter-rotation of the stellar 
and gas components , and the presence of stellar "shell" features in gas-poor systems (e.g. 
Coleman et al.ll2004l ). Diagnosing the occurrence of a previous accretion/merger event in dirr 
galaxies, however, is complicated by their low gravitational potentials, large gas fractions, 
and slow rotation. With slower rotation and low gravitational potentials, their velocity fields 
are more easily disturbed by the energy injection from star formation regions, for example. 

As part of a project t o det ermine what regulates star formation in small galaxies 
(JHunter &: Elmegreenl l2004j . |2006| ). we have obtained multiwavelength data {UBV and Ha 
images) for 136 galaxies, with additional data (e.g. JHK, NUV, FUV) for a subset. These 
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data have recently been augmented by deep, high resolution interferometric atomic hydrogen 
spectral line (H I) observations of 40 of these galaxies as part of the LITTLE THINGS (LT) 
project: Local Irregulars That Trace Luminosity Extremes: The H I Nearby Galaxy Survey 
( |http://www.low ell.edu/users/dah/littlethings/'). As part of our investigation into star for- 
mation in dwarf galaxies, we are interested in looking at "what goes wrong" in BCDs — why 
they have such strong central star formation occurring. VII Zw 403 is one of a few BCD 
galaxies in our sample, and we use it here as a case study for actively star-forming dwarfs. 
Because it is apparently relatively isolated, with no spiral density waves or shear, it is a good 
testbed for examining how star formation might be triggered in active dwarf galaxies. We 
present the results of an initial study here; future papers associated with the LT project will 
present more in-depth analysis and modeling, particularly of the kinematics. 

Although VII Zw 403 has generally been assumed to be a member of the M81 group, 
Karachentsev et al.l (120021 ) state that VII Zw 403, with a distance of 4.4 Mpc, is behind 
the M81 group and is being accelerated towards both it and the Local Group. They find 
that VII Zw 403 is 1. 14 Mpc frorn M81 a nd t hat the M81 group has a zero- velocity radius 
of 1.05 ± 0.07 Mpc. iLvnds et all Jl998h and ISchulte-Ladbeck et all Jl999h determined a 
distance of 4.5 Mpc using the tip of the red giant branch from color-magnitude diagrams 
based on Hubble Space Tele s cope (HST) data. This places it on the far side of the M81 
group. iHunter fc Elmegreenl J2OO4I ) used the NASA/IPAC Extragalactic Database (NED)Q 
to search for nearby galaxies in the plane of the sky within 1 Mpc and ±150 km s~^. The 
only object found (an Im galaxy, KDG 073) had a projected distance of 900 kpc with a 
velocity difference of —32 km s~^. 



Loose fc ThuanI (119851 ) classified VII Zw 403 as an iE galaxy: one with irregular bright 



star formation regions near, but not at the center of, a low-surface brightness envelope. 
Subsequent observations with HST confirmed the presence of a smooth elliptical distribution 



of old, metal-poor red giant stars w ith an exponential fall-off in surface density (JLynds et al. 



1998; Schulte-Ladbeck et al. 



1999 



and referen ces therein). The scale-lengt h for this evolved 



popula tion is about 1 kpc (JLynds et al.lll998l ). which is normal for BCDs (jPapaderos et al. 



1996al ) . There are several bright H II regions with young blue stars loosely cl ustered nearby, 
and diffuse Ha emission as well. Fabry-Perot Ha maps (JThuan et al.l 119871 ) detected little 
rotation, finding only a chaotic velocity field and a velocity dispersion of approximately 30 
km s~^ in the diffuse emission, consistent with stellar winds. 



Initial H I inteferometric observations by iThuan et al.l (120041 ) revealed that although 
the H I is centrally concentrated in general in the galaxy, there are two areas of higher 



^http://nedw"W"w. ipac.caltech.edu/ 
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flux density. The velocity dispersion map showed high values for the region separating 
the two peaks, and they interpret the high dispersion in the region between the peaks as 
due to a lack of flux. As most H I observations of BCDs show smoothly distributed gas 
with only one mostly centralized density peak, the presence of two peaks was intriguing, 
especially if truly separated by a region of lower density. To examine the gas morphology 
and kinematics in greater detail, and to search for possible companions and/or evidence of 
gravitational interactions, we subsequently obtained high spatial and velocity resolution H I 
interferometric data of VII Zw 403 to combine with multiwavelength optical, UV, and IR 
data. These data are presented below. 



2. The Observations 

We have collected images of VII Zw 403 at a variety of wavelengths: Ha, GALEX FUV 
and NUV, ground-based UBVJ, and Spitzer 3.6 /im. The images are shown in Figure [Hand 
observations and data reduction are described below. 



2.1. Optical Images 

We obtained UBV and Ha images of VII Zw 403 with the Perkins 1.8 m telescope at 
Lowell Observatory. The UBV images were obtained on photometric nights with a SITe 
2048x2048 CCD binned 2 x 2 for a field of view of 10.4'. The V and B-band images were 
obtained on UT 13 March 1999 and the ?7-band image on 10 April 1999. The exposure times 
were 3 x 900 s for y, 4 x 1800 s for B, and 3 x 1800 s for U . The binned pixel scale was 
0'.'608, and the seeing on individual and in the final averaged images was 2'.'1 in V , 2'.'8 in 5, 
and 3'.'5 in U . 

The data were reduced using common routines in Image Reduction and Analysis Facility 
(IRAF). Using CCDPROC, we subtracted the electronic pedestal (sometimes referred to as 
the "bias") using the overscan strip, and flat-fielded the images using twilight sky flats. 
Twilight fiats were found, empirically, to work better than do me fiats in removing large- 



scale structure in the CCD response. We used iLandoltl (119921 ) standard stars in 18 fields 
over the 3-night observing run in March 1999 to calibrate the BV photometry and in 13 
fields over the 2-night observing run in April 1999 to calibrate the f7-band photometry. The 
rms of the color terms were 0.015 mag for 5, 0.009 mag for V , and 0.011 mag for U . The 
rms of the zero point terms were 0.015 mag for 5, 0.088 mag for V , and 0.002 mag for U . 
The rms of the extinction coefficients were 0.045 mag for 5, 0.053 mag for V , and 0.005 mag 
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for U. The images in each filter were averaged (IMCOMBINE) in such a way as to ehminate 
cosmic rays while still preserving the photometric integrity of the image. Specifically, images 
were compared pixel by pixel for values that deviated from the average by more than 3 
sigma, based on the CCD noise characteristics. Deviant pixels were then eliminated from 
the average. 

To prepare the images for surface photometry, we edited (IMEDIT) foreground stars and 
background galaxies from the final UBV images. The editing involved replacing a circular 
or rectangular area with an interpolation across the region based on the surroundings. The 
l^-band image was done by hand, and then the cursor log file produced from that was used 
to remove objects from the other two filters. In that way, stars were removed in the same 
way for each filter. We then made a two-dimensional fit to the background (IMSURFIT) 
using the non-galaxy part of the image as seen in a deep display of the image. We subtracted 
(IMARITH) this sky image from the object image to produce a sky-subtracted final image 
for analysis. 

We obtained Ha images of VII Zw 403 on 4 February 1995 using a TI 800 x 800 CCD 
provided to Lowell Observatory by the National Science Foundation. We used the Ohio 
State University Fabry-Perot as a simple 3:1 focal reducer since we were not interested in 
narrow-band spectral information and the Ha filter had a higher efficiency than the Fabry- 
Perot etalon. The Ha filter has a central wavelength of 6566 A and a FWHM of 32 A. The 
contribution from stellar continuum was determined using an off-band filter centered at 6440 
A with a FWHM of 95 A. The Ha observations consisted of two 1300 s images, preceded and 
followed by a 600 s off-band observation. Given the FWHM of the filters, this approximately 
matched the S/N of the on and off-band filters. The pixel scale was 0'.'49, and the seeing on 
the individual and final image is 1'.'7. The electronic pedestal ( "bias" ) was subtracted using 
the overscan strip, and the images were fiat-fielded using dome fiats. 

We shifted, scaled, and subtracted the off-band image from the Ha image to remove 
the stellar continuum and leave only Ha nebular e mission. We calibrated the Ha flux using 



observations of 5 spectro photometric staridard stars (IOkelll974l : IStonelll977l ) and observations 



of the NGC 2363 nebula (JKennicutt et al.lll980l ) taken over the course of the 6-night observing 



run. These two types of calibrations agreed with each other to within 4%. The Ha flux has 
been corrected for the shift of the bandpass with temperature. The central wavelength shifts 
blueward by 0.018 nm for each degree of temperature decrease from 20° C, and so this needs 
to be considered in calibrating the flux through these narrow-band fllters as it sometimes 
moves the object emission line off the peak of the transmission function. The Ha flux has 
also been corrected for contamination by the [NII]A6584 emission line in the bandpass (a 1% 
correction) . 



2.2. Near-IR Image 

We obtained the J-band image of VII Zw 403 with the Perkins 1.8 m telescope at Lowel l 



Observatory with the Ohio State Infrared Imager- Spectrograph (OSIRIS iDePoy et al.lll993[ ) 
on the photometric night of UT 2 June 1996. The instrument consisted of four mosaiced 
detectors to create a total array of 256x256 pixels with a field of view of 6.4'. The telescope 
was nodded to place VII Zw 403 first in one half of the array and then in the other half, 
so that sky could also be measured in the half without the galaxy. The galaxy was offset a 
few arcseconds between images in order to allow removal of hot pixels when the frames were 
averaged later. The image of VII Zw 403 is a combination of 48 exposures of 26 s each. The 
pixel scale was 1'.'50, and the seeing of the final image was 3'.'2. 

The data were reduced in IRAF. Zero second bias frames were subtracted from all of 
the data, the image pixel values were corrected for non-linearity effects, and pixel to pixel 
sensitivity variations were removed using observations of a white screen in the dome. The 
non-linearity corrections were determined using a series of dome flats of increasing exposure 
time, to 29 seconds in steps of 1 second, sandwiched between darks. 

Eight UKIRT standard stars from the list compiled by S. Courteau (1996, private com- 
munication) were used to calibrate the photometry. Each star was observed 4 times with 
the star placed in each quadrant in turn, and sky for a specific observation was taken from 
the observations with the star in the other quadrants. The rms of the zero point of the 
J-band calibration was 0.12 mag and of the extinction term was 0.08 mag. The California 
Institute of Technology (GIT) photometric system is 8% bluer in J — if and 4% bluer in 
H — K compared to the UKIRT system. The star for ming region of VII Z w 403 shows 



up faintly on the Two Micron All Sky Survey (2MASS ISkrutskie et al.l |2006| ) J-band im- 
age: aJ_asky_990319n0730056.fits. On the 2MASS image, the flux in a 22" radius aperture 
centered on the star-forming regions has a magnitude of 13.73±0.06. Our J-band image 
gives 13.71±0.01. In spite of the difference in calibration systems (2MASS defines their own 
system), the magnitudes are the same within the uncertainties. 



2.3. Ultraviolet Data 

We obtained near-UV (NUV) and far-UV (FUV) images of VII Zw 403 from the Galaxy 
Evolution Explorer (GALEX) archivecl. The images were part of the Nearby Galaxy Survey. 



2 



Based on observations made with the NASA Galaxy Evolution Explorer (JMartin et ahl 120051) . GALEX 



is operated for NASA by the Cahfornia Institute of Technology under NASA contract NAS5-98034. 
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The FUV filter has a bandpass of 1350-1750 A, an effective wavelength of 1516 A, and a 
resolution of 4'.'0. The FUV exposure time was 1617 s. The NUV filter has a bandpass of 
1750-2800 A, an effective wavelength of 2267 A, and a resolution of 5'.'6. The NUV integration 
time was 1479 s. The images were processed through the GALEX GR2/3 pipeline and were 
retrieved as final intensity maps with a 1'.'5 pixel scale. The GALEX field of view is a 
circle with 1?2 diameter, and we extracted a portion around VII Zw 403. We subtracted 
(IMEDIT) or masked foreground stars and background galaxies, and subtracted a low order 
two-dimensional fit to the background (IMSURFIT). The decision of what objects did not 
belong to the galaxy and should be removed was based on morphology (extended or stellar 
profiles), color, and similarity with other objects in the field well beyond the galaxy itself. 



2.4. 3.6 /im Data 

For completeness and to facilitate comparison, we include the 3 .6 /im image taken of 



VII Zw 403 with the InfraRed Array Camera (IRAC lFazio et al.lll982l ) on the Spitzer Space 



Telescope {Spitzer). The data were obtained from the Spitzer archive and red uced with the 



Basic Calibrated Data pipeline. The full treatment of the data is described by lHunter et al. 



(120061) 



2.5. H I Data 

21 cm spectral line observations of the galaxy were made using the B, C, and D con- 
figurations of the NRA0l2|. Very Large Array ( VL A) . The B and D-configuration data were 
taken by Deidre Hunter; the C configuration by Trinh X. Thuan (proposal code AH0453), 
who kindly supplied us with these observations. The D configuration, with its short spac- 
ings, provides high sensitivity but a relatively low spatial resolution of ~ 1'. The C and B 
configurations, with longer baselines, provide higher resolution (nominally 12'.'5 and 4", re- 
spectively) but lower surface brightness sensitivity. After editing and calibration (discussed 
below), the three data sets were combined to produce images with both high sensitivity and 
high angular resolution. 

Time on source for the C and D configuration observations was approximately 4 hours 
each; for the B, approximately 9 hours. Observations for all configurations were made using 



■^The National Radio Astrononiy Observatory is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 
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a correlator setting using 2IF mode with Hanning smoothing with a total bandwidth of 1.56 
MHz. This results in 128 channels with a channel separation of 12.2 kHz (2.6 km s~^). For 
spectral line data that are taken using on-line Hanning smoothing, the velocity resolution 
is equal to the channel separation, and so is also 2.6 km s~^. Observational information is 
provided in Table [H 

The D configuration observations were set up to be consistent with the previously ex- 
isting C configuration observations, so the two sets of observations could be combined after 
calibration for analysis. The central observing frequencies and velocities were the same, other 
than the usual changes in sky frequency due to the day, time, and location of the sources on 
the sky. However, the B configuration observations were taken during the beginning of the 
transition of the VLA to the Extended VLA (EVLA); so doppler tracking while observing 
was not available. The observations were done at a fixed frequency which corresponded with 
the central velocity (—60 km s~^) from the C and D observations to within 3 or 4 channels 
(approximately 12 kHz). Before the B data were combined with the other two data sets (see 
below), the AIPS task CVEL was used to shift the spectra in the B data set to correct for the 
Earth's rotation and motion. Careful comparison of the resulting spectrum after applying 
CVEL with those from the C and D data ensured that all three data sets were now consistent 
in frequency and velocity space. 



2.5.1. Calibration 

Calibration of all three data sets was performed independently, using the standard rou- 
tines in AIPS. For the D configuration data, the flux and bandpass calibrator was 1411-1-522 
(3C 295; all sources are in J2000) and the phase calibrator was 1435-1-760 (J2000). The C 
configuration observations used 1044-1-809 for phase calibration and 1331-1-305 (3C 286) for 
the flux and bandpass calibration. The phase calibrator for the B configuration data was 
the same as for the D (1435-1-760), and the flux and bandpass calibrator was 0137-f331 (3C 
48). 

The data reduction process was complicated by the presence of both absorption and 
emission from the Milky Way as VII Zw 403's systemic velocity is ^ —100 km s~^ . The 



single-dish spectrum (Figure |2]) from the Greenbank 91m telescope ( jThuan fc Martinlll98l[ ) 
shows a broad absorption feature from ~ —55 to —30 km s~^, with a strong emission feature 
from —30 to ~ -|-5 km s~^. That emission feature shows a dip due to a narrow absorption 
line at a velocity of about km s~^. 

To reduce the amount of data potentially contaminated by Galactic H I in the D configu- 
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ration calibrator data, a pseudo-continuum "channel 0" dataset was constructed by averaging 
together the inner 50% of the line data (normally, the inner 75% is used). This probably 
did not remove all the emission from the Milky Way, but was the best compromise given the 
situation. Because these were daytime observations, solar contamination was also present on 
the shortest baselines for the calibrators as well as for the source data. The affected short 
baselines for the calibrator data were not used during the calibration by the application of an 
appropriately chosen "uvrange" (only baselines longer than 0.55 klambda were used). Once 
a satisfactory calibration was achieved, it was applied to the source line data. The band- 
pass calibrator was only affected by Milky Way H I in a few channels, and as those hadn't 
stood out as unusual in the previous solutions, the bandpass calibration was done, again 
applying the uvrange limit (0.55 kA) for the solar contamination. The continuum emission 
was removed using UVLIN, which does a linear fit to specified (line-free) channels and then 
subtracts the fitted baseline from the spectral line channels. The channels used for the fit 
were chosen to be free of both line emission and Galactic contamination. 

For the C configuration observations, the continuum source 1331+305 (J2000) was ob- 
served for use as a flux and bandpass calibrator; the nearby continuum source 1044+809 
(J2000) was observed for phase calibration. There was little Galactic or solar contamination 
present in these data. For the C data, calibration proceeded normally. During editing of the 
channel source data, amplitudes above 0.4 Jy were clipped to remove anomalously high 
amplitude visibilities on the shortest baselines during the beginning of the observation run. 
The editing and calibration tables were copied over to the line data, and a dirty map of 
the source data was made to check for interference and look for line-free channels to use for 
continuum subtraction. Continuum subtraction using the line and inteference-free channels 
thus identified was done in the u-v plane using the task UVLSF, which works similarly to 
UVLIN. 

The B configuration observations used 1331+305 (J2000) for flux and bandpass calibra- 
tion, and 1435+760 (J2000) for phase calibration. The data from four antennas were flagged 
as unusable during the data collection process due to EVLA testing; four more antennas 
were not in use during the observations. During calibration, several large phase jumps in the 
phase calibrator data, on the order of 100 degrees, were observed for some of the antennas. 
All of these jumps occurred between scans near the end of the observing run. Changing 
the reference antenna for the calibration did not change the number or magnitude of the 
problematic jumps. After further consideration and testing, the galaxy source data from the 
scans that were in between the calibrator scans with phase jumps larger than 40 degrees 
were removed by flagging. Careful examination of the resulting images revealed no artifacts 
that would indicate phase errors with the data, such as asymmetric stripes, rings, or other 
odd-looking features. Self-calibration was also tried, but the resulting images were essen- 
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tially identical to those from the edited data. Resuhs here are from the edited data that 
were not self-cahbrated. As with the D data, continuum subtraction was done in the u-v 
plane using UVLIN. There was no apparent contamination from Milky Way emission visible 
in the B configuration data. 

To check the calibration and data quality, the data from all three configurations were 
imaged separately using the task IMAGR, which uses a fast-fourier transform to create a cube 
of images, one image for each frequency (velocity) channel. During imaging, the effects of 
side-lobes of the non-gaussian features of the "dirty" beam were removed using the clean 
algorithm (with a flux cut-off equal to the r.m.s. in a line- free channel). With the exception 
of the Galactic contamination obvious in the D configuration data (and faintly in a few 
channels of the C data), all three data sets were deemed acceptable. 

A spectrum from our D configuration data (uniform weighting, with a robust factor 
of +l)u is shown in Figure |3l Some of the Galactic contamination is strongly visible around 
km s~^, which doesn't overlap with VII Zw 403's emissio n, visible around —100 km s~^. 



Many of the same features are present here as are in the iThuan fc Martini (1198 ll ) single- 
dish spectrum, although at different levels. Presumably this is due to the spatial scale of 
the features; the D configuration is "blind" to structures larger than about 15' per spectral 
channel. 

The D configuration data were then reimaged using natural weighting to emphasize faint 
emission; the spatial resolution is 72'.'9 x 53'.'6 and the la noise in a single channel is 0.93 
mJy beam~^ (see Tabled]). Our integrated flux from the D configuration data is 11.5 Jy km 
s~^, which is 25% lower than previous single-dish measurements. (We did not correct for the 
sensitivity fall-off of the primary beam, as the maximum extent of H I emission was approx- 
imately 6', which is well within the 30' full-width at half-maximum ( FWHM) of the primary 



beam at 21 cm.) Estimates from s ingle dish observation s are 14.6 (iHuchtmeier fc Richter 



1986h . 15.8 jThuan fc Martini 1 1 98 ih . and 16.2 Jy km s'^ Jfluchtmeier fc Richteil EJ ) . We 



checked our flux calibration by comparing the fluxes of continuum so urces from our data to 
the catalogued fluxes listed in the NRAO VIA Sky Survey (NVsfl; ICondon et all Jl998h l 



There are two catalogued continuum sources close (3.7' and 3.5') to the galaxy; our measured 
fluxes are 79% and 76% of the NVSS fluxes, respectively. If we adjust our D configuration 
integrated flux by 20%, this still only comes to 13.8 Jy km s~^. A comparison of a spectrum 



''The ROBUST parameter allows the user to choose a weighting scheme between fully uniform (—5) and 
fully natural (+5). This produces images with either higher spatial resolution but poorer signal-to- noise 
(fully uniform) or lower spatial resolution but higher signal-to-noise (fully natural). 

^http://www. cv.nrao.edu/nvss/ 
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through the D configuration data cube restricted in area to the spatial ext ent of the galaxy 



on the sky is shown in Figure HI along with the single-dish spectrum from iThuan &: Martin 



( ll98ll ). Presumably, all the single-dish integrated fluxes include the second emission peak 
located at about —70 km s~^. We have determined that this secondary emission peak is 
almost certainly from the Milky Way, not VII Zw 403 itself. Comparing just the emission 
from VII Zw 403, it is clear that our observations have detected all of the flux seen in the 
single-dish spectrum. 



This is a different conclusion than that reached by iThuan et al.l ( 120041 ) . They present 
the VLA C configuration data that we show here, but did not have the D configuration 
data and so were not aware of the Galactic emission in the field of view. They also present 



a comparison of their C data spectrum to the original iThuan Sz Martini (Il98l[ ) single-dish 
spectrum. Their C spectrum resembles our D spectrum without the Milky Way emission; 
indeed their total integrated flux from the C configuration data is 11.0 Jy km s~^, which 
is close to our 11.5 Jy km s~^ value. However, without knowing that the single-dish flux 
included the Galactic emission as well, they were forced to conclude that their VLA data 
were missing some (presumably large-scale) emission from VII Zw 403. 



2.5.2. Imaging 

Channel maps from part of the D configuration data cube are shown in Figure |5l The 
Galactic contamination is clearly visible in the channels from —23.9 to —78.0 km s~^. Several 
of these channels also contain emission from VII Zw 403. Before integrating the cube to 
produce the moment maps showing the integrated flux, velocity field, and velocity dispersion, 
the affected channels were selectively blanked to remove the Galactic emission, leaving only 
the small region in the center of each field that contained emission from VII Zw 403. This 
partially blanked cube was then integrated in velocity using the channels from —60.0 to 
— 147.5 km s~^ with the AlPS task momnt. To avoid integrating up noise as well as signal, 
MOMNT creates a mask cube smoothed in space (by convolving with a gaussian twice the 
size of the beam) and velocity (by banning smoothing); pixels falling below the specified flux 
cutoff in the mask cube are then blanked in the original cube prior to integration. For this 
dataset, a flux cutoff of twice the single-channel r.m.s. was used. These moment maps were 
used for the analysis discussed in § HI 

To produce a data set with good sensitivity and enhanced but still moderate resolution, 
the data from the C and D configurations were combined in the u-v plane using the AlPS 
task DBCON and then imaged using the AlPS task IMAGR with a ROBUST weighting factor of 
+1 (we will refer to this cube as the C+Di cube). This weighting was chosen based on tests 
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to determine the best balance of resolution and flux recovery (sensitivity) for the purposes 
of velocity field analysis, primarily. Although the Galactic contamination wasn't as strongly 
present as in the D-configuration-only data, it was still present at high enough levels to 
warrant concern. We therefore blanked everything but the emission from VII Zw 403 (in the 
center of the field) in the eight affected channels before integrating. A flux cutoff of twice 
the single channel r.m.s (see Table [2]) was used to produce the combined data moment maps. 
This data set was used for the velocity analysis and H I surface density analyses, below. 

The data from all three configurations were also combined in the u-v plane and imaged. 
Galactic contamination was not detectable in any of the resulting data cubes, so no extra 
blanking was necessary. We tested a variety of ROBUST weighting factors while imaging to 
determine the best values to provide the desired resolutions while recovering the most flux 
possible for each resolution. After extensive testing and comparison, we created three data 
cubes with varying sensitivities and resolutions: a robust of produced a medium resolution 
data set (B+C+Dq) which was used to examine the morphology of the integrated flux (the 
first moment map); a ROBUST of +1 produced a lower resolution data set (B+C+Di) used 
to analyze the velocity dispersion in the gas (the second moment map); and a ROBUST of —1 
was used to create a high resolution (B+C+D_i) image of the integrated flux in the highest 
density regions of the galaxy. Details (resolutions, noise levels) are given in Table |2j 



3. The Optical View of the Galaxy 

The images of VII Zw 403 are shown in Figure [H The galaxy is elliptical-shaped with a 
moderate degree of ellipticity implying a relatively high inclination angle. At outer l^-band 
isophote levels VII Zw 403 is quite symmetrical. However, at intermediate isophotes (radius 
40-60") there is an extension of stars to the north that is not mirrored to the south. This 



type of feature is seen in other irregulars as well, such as NGC 2366 ( JHunter et al.ll200ll ) 



3.1. Broad-band Surface Photometry 

From the V^-band image, we determined the morphological center of the galaxy, b/a, 
and position angle of the major axis from a contour plot of the image block-averaged by 5 to 
increase the signal-to-noise in the outer parts. We used an outer contour corresponding to 
a nv, not corrected for reddening, of 27.2 mag arcsec"^. The major axis and position angle 
were determined to be the long line of symmetry around which the galaxy could be folded. 
The center was the mid-point bisecting the isophote along the major and minor axes. 
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The center of the galaxy is at 11*^ 24"^ 33'.'4, 79° 16' 10" (B1950), with an uncertainty 
of 2" in eac h direction. We also determined the center of the galaxy from Spitzer IRAC 3.6 
/^m images (JHunter et al.ll2006l ). with an estimated uncertainty of 3". The 3.6 /im center is 
offset from that of the l^-band by 4" W and 3" N, and this then more realistically describes 
the uncertainty in the position of the center of the stellar body. Below we will derive the 
H I kinematic center as well. That center is at 11*^ 24" 35n±10!3, 79° 16' 18±27" (B1950). 
The H I kinematic center is offset from that of the V^-band by 5" E and 8" N, well within the 
large error bars of the H I kinematic center. For the surface photometry we used the V^-band 
morphological center. We measured the minor-to-major axis ratio b/a on the l^-band image 
to be 0.49. Dwarf irregular galaxie s are considered to be th i cker than spiral disks and an 
intrinsic (b/a)o o f 0.3 is suggested (JHodge &: Hitchcock! Il966l : Ivan den Berghlll988l : but see 
Sung et al.lll998l ). This b/a then implies an inclination of 66° for VII Zw 403. We measured 
the position angle of the major axis to be —11°. 

The surface photometry (ELLIPSE in STSDAS) used 10 ellipses starting at a semi-major 
axis of 9" and progressing to 91" in steps of 9'.'1. All of the other images were geometrically 
transformed (GEOMAP, GEOTRAN) to match the l^-band image and elliptical surface 
photometry was done identically on all images. The geometric transformation used stars 
identified in both the l^-band image and the target image to determine the surface distortions 
necessary to make the target image have the same size, pixel scale, and orientation of the 
l^-band image. 

The surface photometry is collected in Figure O The surface photometry and colors 
have been corrected for reddening using a total E(B — V)t=E(B — V^)/-|-0.05, where the 



foreground reddening E{B — V)f is 0.023 (JBurstein fc Heileslll984l ). An internal reddening 
of 0.05 in E{B — V) was used to be consistent wit h Balmer decrement measurements of H II 
regions in dim galaxies. iHunter fc Hoffman! (119991 ) found an average E{B — V) of 0.1 for H II 
regions in 39 dim galaxies. The stars outside of H II regions were taken to have a reddening 
of half of this value. We used Burstein and Heiles values for the foreground reddening rather 
than those of !Schlegel et al.! (11998! ) because the latter values are systematically higher than 
the former and some are inconsistent with Balmer decrement measurements. Our value of 



the reddening of the stars in VII Zw 403 is close to that of iLynds et al.l (11998! ) who obtained 
an E{B — V) of 0.04-0.08 for stars and stellar associations f rom analysis of color-magnitude 
diagrams. We use the reddening law of ICardelli et al.l (11989! ) and Av/E{B — V)=3.1 for the 
optical, near-IR, and FUV reddening corrections. The NUV reddeni ng is affected by th e 
2175 A bump. For the NUV, we have adopted the reddening formula of IWyder et al.l (120071 ). 



We determined R25, the radius at a i?-band surface brightness of 25 mag arcsec ^, 
to be 36'.'5 which is 0.78 kpc at the galaxy. Our radius is 10% smaller than that given by 
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de Vaucouleurs et al.l (Il99l[ ). The Holmberg radius, Rn, originally defined to a photographic 
surface brightness, is measured at an equivalent 5-band surface brightness fiB = 26.7 — 
0.149(5 — V). For our measured {B — V)q of 0.3, the Holmberg radius is determined at 
a hbq of 26.66 magnitudes arcsec"^. We find that R^ is 69" which is 1.47 kpc at the 
galaxy. Both R25 and Rh are determined from the reddening corrected surface photometry. 
Integrated properties of VII Zw 403 are collected in Table [31 

We fit an exponential disk to the broad-band photometry beyond 0.2 kpc, and the fits 
to the V^-band, 3.6 /im, and NUV data are shown in the upper panel of Figure El The 
outer disk of VII Zw 403 is fit well with an exponential disk profile having a central surface 
brightness /i^ of 22.94±0.09 magnitudes arcsec"^ and a disk scale length R^) of 0.50±0.01 
kpc. Interior to a radius of 1.0 kpc, the surface brightness profile is steeper, although it too 
is fit well with an exponential disk with a disk scale length of 0.37±0.01 kpc. The 3.6 /im 
surface brightness profile follows the slope of the V^-band reasonably well while the NUV 
profile is considerably steeper. The similarity between the V and the 3.6 /im profiles implies 
that the star formation as a function of radius integrated over a 1-Gyr time scale is the same 
as that integrated over the lifetime of the galaxy. The steeper NUV profile suggests that 
star formation over the last 200 Myrs has been concentrated to the center of the galaxy. 

Surface brightness profiles of BCDs are generally found to contain multiple components 
with the outer component representing the galaxy underlying the inte nse current star forma- 
tion (jPapaderos et al.lll996bt ICairos et al.l 120031 : iNoeske et al.ll2003[ ). The central portions 
of these galaxies, on the other hand, are dominate d by the current s tar for mation, and the 
profiles are often steeper or more complex there. iPapaderos et al.l (jl996al ) found that the 
scale-length of the outer disks of BCDs are sm aller, by half on average, than those of Im 
galaxies. However, JHunter fc Elmegreenl ( l2006l ) find for a sample of 94 Im systems and 24 
BCDs that the scale-lengths of the BCDs are only slightly offset (in the sense of shorter 
scale-length at comparable My) from those of most Im galaxies. It is the higher scale-length 
tail in the number distribution that is more drastically reduced for their BCD sample com- 
pared to dims. Compared to these samples, VII Zw 403's scale-length lies at the peak in the 
distribution of dims and BCDs. 



Integrated colors of VII Zw 403 are given i n Table El In a UBV color - color diagram for 
a large sample of Im, BCD, and Sm galaxies ( jHunter fc Elmegreenll2006l ). VII Zw 403 lies 
at the blue side of the range, clearly dominated by a young stellar population. On the other 
hand, the integrated {V — J)n places VII Zw 403 i n the middle of the range shown by Im 
and BCDs in that color ( jHunter fc Elmegreenll2006l ). 



Azimuthally-averaged radial profiles of the colors are shown in Figure |6l In most Im 
galaxies the azimuthally-averaged colors are generally constant over the disk of the galaxy; 
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irregular galaxies do not usually have color gradients. VII Zw 403, however, does show a 
gradient towards redder colors in the outer part of the galaxy. {B — V)o and (f/ — B)q track 
each other and change by about 0.5 mag from the center into the outer galaxy, beyond the 
Holmberg radius. We can only trace {V — J)q to 1.1 kpc, but it is approximately constant 
over this radial range. {FUV — NUV)o also is roughly constant throughout the optical 
galaxy to R25. This trend in the colors is due to the fact that star formation is concentrated 
to the center of the galaxy. Stellar populations becoming redder as one goes outwards from 



the c enter is seen in many BCDs (jPapaderos et al.lll996bl : ICairos et al.ll2003t iNoeske et al. 



2003|). 



We have modeled the UV and optical colors of VIIZw403's annular surface photometry 
using population synthesis fitting. These mode l s use simple evolutionary hist ories and the 



stellar populations library of lBruzual &: CharlotI (|2003[ ). Details can be found in lHunter et al. 



(J2OIOI ). We find that the colors of the central region (<0.5 kpc) is dominated by a stellar 
population with an age of about 400±200 Myrs, while the colors beyond that radius are fit 
with stellar populations of 1-3 Gyr of age. 



3.2. Star-Forming Regions 

Ha traces star formation over the past 10 Myr, and Figure [T] shows that most of the 
current star formation in VII Zw 403 is found in one large complex located 6" east and 2" 
north of the center of the galaxy, only 130 pc from the nominal center determined from the 
outer l^-band isophote. This type of central concentration of star formation is often seen in 
Blue Compact Dwarfs, and is probably responsible for VII Zw 403's classification as a BCD. 

We show the azimuthally-averaged Ha surface brightness in Figure |6l Although the 
current star formation in VII Zw 403 is contained in a giant H II complex at the galaxy 
center (the Ha emission seen in Figured]), because the H II complex is large. Ha emission 
appears in elliptical annuli out to a radius of 1.3 kpc. The Ha surface brightness drops 
off faster than the starlight in the \^-band and 3.6 fim image, dropping only a little faster 
beyond 0.75 kpc and ending at 1.3 kpc. 

The FUV and NUV passbands measure the light directly from stars and are dominated 
by young massive stars formed over the past 200 Myr. Thus, the UV is sensitive to star 
formation over a longer timescale than Ha. In Figure E] the UV profile, however, follows the 
Ha profile closely. Thus, the ultraviolet surface brightness profile is dominated by the same 
central star-forming event as we see in Ha (see Figured]). 

The integrated Ha and UV luminosities and derived star formation rates of the galaxy 
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are given in Ta ble [3l The star form ation rates derived from Lna and the FUV lumin osity use 
the fo rmulae of iHunter et al.l (120101 ) that integrates from 0.1 Mq to 100 Mq with a lSalpeter 
(119551 ) stellar initial mass function. The star formation rate derived from Ha is similar to 
that derived from the UV. VII Zw 403 has a star formation rate per unit area that is norma l, 
but towards the high end of the range, for irregular galaxies (JHunter fc Elmegreenll2006l ). 
as shown in Figure [71 Again, this is consistent with the galaxy's classificatio n as a BCD. 



Howev er, modeling of the color-magnitude diagram of the stellar population by iLynds et al. 



(119981 ) suggests that the star formation rate was much higher about 600-800 Myr ago than 
it is today, by as much as a factor of 30. 

The V and 3.6 /xm-band radial profiles in Figure [6] are smooth exponential disks with 
about the same slope in the central part of the galaxy, but the Ha and NUV are peaked in 
the center and fall faster than V or 3.6 yum. This indicates a change in the primary site of 
star formation about 200 Myr ago, the timescale the NUV is sensitive to. This is a relatively 
short dynamical time (~ two galactic rotations; see § 14.31) for this slowly rotating galaxy. 
We cannot determine the cause of this shift in star formation location, but it is unlikely to 
have been caused by a collision or accretion event. The timescale is short enough that we 
should still see either two galactic nuclei or tidal debris, which we don't. A central collapse 
to a bar could explain it, but we see no hint of a bar in the V, J, or 3.6 /xm images. 

At its current rate of consumption, the galaxy can turn gas into stars for another 5 Gyr if 



all of the gas associated with the galaxy can be used. The ti mescale to run out o : 



10 Gyr if recycling of gas from dying stars is also considered (JBrinchmann et al 



times cales to exhaust the gas are somewhat short for irregular galaxies (JHunter fc Elmegreen 



gas b ecomes 
20oi). Thes e 



20041 ) , but nevertheless, the galaxy is in no danger of running out of the fuel for star formation 
any time soon. 



4. HI Results 



As H I is the raw material from which molecular clouds, and then stars, form, the H I 
properties of a galaxy are fundamental to understanding its past, current, and future star 
formation. The total H I mass from the integrated D configuration data is 5.23 x 10^ Mq, 
which is below the median of 3.0 x 10*^ for the large sample of star-forming galaxies studied by 
Salzer et al.l (120021 ). This value for the hydrogen mass includes only the emission we believe 
is associated with VII Zw 403 itself. This does not include the secondary emission peak seen 
in the spectrum around —70 km s~^ that we believe is due to Galactic contamination (see 
§[231]). 
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The Mnr/Ln i s 0.91 This is within the range (0.2-2.0) of the sample of 11 BCDs 
from Ivan Zee et aLl (jl998l . 1200 ll ). but is on the high end of the distribution for the 24 
BCDs (~ 0.005 to ~ 1.8) , and below the median for Ims (~ 1.3) shown i n Fig ure 4 of 



Hunter fc Elmegreenl (120041 ) . For comparison with the sample in ISalzer et al.l (120021 ) , we fol- 
low their method to "correct" the global luminosity by removing the effect from the starburst 
so that it is more representative of the bulk o f the stars in the g alaxy. This increases the 
Mhi/L-q ratio to 1.82, which is higher than the ISalzer et al.l (12002! ) median of 1.3. Although 
the H I mass for VII Zw 403 is relatively low on an absolute scale, its relatively high Mm/ L-q 
(even "uncorrected" for the starburst) could indicate that it is either a somewhat gas-rich 
BCD or is currently "luminosity-poor" as it is presumably po st-starburst: it sho uld have 
been brighter during the burst of 600-800 Myr ago detected by lLynds et al.l (119981 ). 



4.1. HI Surface Density 

The azimuthally averaged surface density profiles of the H I gas for VII Zw 403 using 
the center, position angle, and inclination from both the optical and H I data are shown in 
Figure [HI The profiles were calculated by integrating the H I from the C+Di map (beam size 
34'.'1 X 27^6) using 30" rings with (1) a center position of 11^^24'^33!5, 79° 16' 10'.'3 (B1950.0), 
a position angle of —11°, and an inclination of 66° taken from the ^-band photometry (§IS]); 
and (2) a center position of 11^24" 35!1 ± 10!3, 79° 16' 18 ± 27"(B1950.0), a position angle 
of 47°, and an inclination of 77° taken from the NE half of the H I rotation curve analysis 

(SD- 

Several Im galaxies are also plotted in Figure |H1 The galaxies were chosen to have 
similar linear resolutions in that the ratio of the beam size to either the disk scale length 
{Rt)) or the half-light radius (-R1/2) is approximately the same. The ratios range from 1.1 
to 2.4 for beam/i?D and 0.5 to 1 for beam/i?i/2- The values for VII Zw 403 are 1.4 and 1.2, 
respectively. We can see from the figure that the H I surface density for VII Zw 403 derived 
using the H I parameters falls off smoothly and slightly more steeply in the central regions 
than those for the Ims; more so when using the optical parameters. I n corn parison to the 
BCD and dirr H I surface density profiles presented in Ivan Zee et al.l ( 1l998l ). the shape of 
VII Zw 403's profile is similar to those of the BCDs, and somewhat steeper than those for 
the dirrs. The central surface densities from both the optical and H I parameters for VII Zw 
403 are seen to be low compared to the BCDs in the sample; when using the H I parameters, 
the central surface density is like that of the dIrrs. 
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4.2. Morphology 

Figure [9] shows the integrated C+Di H I contours on the V image, with the contrast 
set to show the full extent of the galaxy's optical emission. This figure shows that there is a 
large, low-level irregular skirt of H I gas that extends beyond the optical emission, especially 
toward the south. Irregular gas distributions that extend b eyond the optical emission, which 



are common in field galaxies, a re also found in other BCDs ( jTaylor et al.lll994J : Ivan Zee et al. 



19981 . I2OOII : iThuan et al.ll2004j ). From the D configuration data (not shown), the total extent 



(diameter) of the H I {2Rhi) out to the 1 x 10^^ c m ^ lev e l is 5 .24 (deconvolved). This is 



2.3 -Rh; which is typical for the Ims in Figure 13 of iHunten (119971 ). 



At the spatial resolution of the C+Di image (~ 700 x 600 pc), the H I appears smoothly 
distributed with the density increasing t oward the center of the galaxy. This smoothly 



increasing density i s often seen in BCDs ( iTaylor et al.l Il994j : ISimpson &: GottesmanI 12000 



van Zee et al. 



1998 



200 if): and although some "dumpiness" of the H I is also seen in 



BCDs ( Ivan Zee et al.lll998l ). it is not like the highly-structured gas di stribution often ob- 



served in dirrs, where shells, fi 



Puche et al. 



1992 



Stewart et al. 



ilYoung fc Lo 



2000; 



Qtt et al. 



aments, and loops are frequently visible (iMeurer et al. 



19961 : Martini 



1998 



2001t ISimpson et al. 



Meurer et al. 1998; Walter fc Brinks 



19921; 



199ft; 



2005aU bl) . Comparison of the H I distri 



but ion in VII Zw 403 to those in dIrrs with similar linear spa tial resolution, such as DDO 
88 ( ISimpson et al.ll2005al ) and DDO 43 (ISimpson et al.ll2005bl ). shows that at these scales, 
this increased "smoothness" in VII Zw 403 is real. Clumping of the H I on the scales seen 
in dIrrs is not visible in VII Zw 403. 

At higher spatial resolution, however, some structure begins to be revealed in the H I. 
Figure Uni shows the integrated H I fiux from the B+C+Dq data cube, with a linear resolution 
of 200 xl60 pc (9" X 8"), as compared to the 750 x 600 pc resolution of the C+Di image. 
(Please note the change in field of view of the two images; the B+C+Dq image shows only 
the inner 2', which corresponds roughly to the area covered by the central three contours 
in the C-I-Di map shown in Figure [HI) This image shows that the central H I distribution 
contains a region of lower density separating a large clump to the south of the center and a 
much smaller clump to the north, connected by a ridge on the west side. This low density 
region resembles a partial cavity or hole, with decreasing density to the east side. This is 
further discussed below. 
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4.3. The Velocity Field 



In Figure [TT] we show the channel maps from the C+Di data cube that were integrated 
to produce the flux-weighted H I velocity fleld shown in Figure [121 and in Figure [13] as 
isovelocity contours superposed on the V-band image of VII Zw 403. One can see here that 
there is ordered rotation, with the NE portion of the galaxy approaching and the SW half 
of the galaxy receding. However, it is also immediately apparent that the velocity fleld is 
complex. First, there is a curve in the isovelocity contours so that the NE and SW portions 
of the galaxy have distinctly different kinematic axes. The kinematic axis of the NE half of 
the galaxy is aligned roughly at an angle of +45°. The SW part of the galaxy, on the other 
hand, appears to have a line of nodes at about +25°, which is aligned more closely with the 
optical galaxy's major axis at —11°. Second, there are places where the velocity fleld is not 
well behaved. In particular on the east side of the galaxy the isovelocity contours curve back 
upward in a complex manner, and in the south there are "islands" of different velocities. 

As strange as this velocity fleld appears, there are other dwarf galaxies with curves in the 
position angles of their isovelocity contours similar to what we see in VII Zw 403. A review of 



the li terature re vealed ten other BCDs wi th this type of featur e (ESO338-IG04: lOstlin et al. 



2001; Haro 2: 



NGC 1569: Stil fc Israel 



Bravo-A 



aro et al.l I2004J : Ha ro 33, Mrk 51: ISimpson fc Gottesman 



2002b; NGC 3738: Stil fc Israel 2002a; NGC 1705: 



2000 



Meurer et al. 



19981 : UGC 4483, UM462: Ivan Zee et al.l [19981 UM323: Ivan Zee et al.lboOli ). and there are 



undoubtedly others. However, this phenomenon is not exclusively that of BCDs. We also 
found a coraparabl e number of Im galaxies wit h this sort of feature (DDO 22, DDO 68: 
Stil fc Israeli l2002al: DDO 69: IVoung fc Lolll996l : DDO 165, DDO 187, UG C 3698: ISwaters 



Cote et al.ll2000l : NGC 2366 



Hunter et al 



I999I: DDO 2 10: iBegum fc Chengaluil booj: DD O 216: IVoung et al.ll2003[ ESO 245 -G005 



200lt NGC 4163: ISimpson fc Gottesmanll2000h . 



Since some of the Im galaxies have no current star formation activity or very little, we cannot 
associate this velocity fleld feature with the starburst that makes VII Zw 403 special. In 
addition, this is not relegated only to slow rotators, such as VII Zw 403. In the sample 
taken from the literature, maximum rotation speeds range from very small up to 69 km 
s~^. Distortions in velocity flelds are often blamed on gravitational interactions with other 
galaxies. In these galaxies, some are known to have small companions and others are not. 

To get some idea of the rotation of VII Zw 403, using the task gal in AlPS we flt the two 
halves of the velocity fleld separately and began with the NE or approaching half. We used 
the C+Di moment map which has a beam size of 34'.'lx27'.'6. With the H I peak a s a flrst 



guess on the position of the kinematic center, we flt the fleld with a Brandt function ([Brandt 



1960[). We also experimented with changing the initial conditions: we varied the position 
angle from 34-65°, the inclination from 66°to 72°, and maximum rotational velocity from 13 
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tp 15 km s^^. We also varied the radial extent of the fit to raidd of 130", 90", and 75"; the 
center up to eight pixels; and the systemic velocity by a few km s~^. We tried fitting the 
whole field and just one side at a time while varying initial parameters. For the final rotation 
curve (s) we used a position angle of 40°, an inclination of 66°, and a rotational velocity of 
13 km s~^ for the initial guesses. From this we determined the center (ll'^ 24™ 35!1±10'^3, 
79° 16' 18±27", B1950.0) and systemic velocity ( — 103 ± 1 km s~^). The experiments with 
varying initial guesses suggested that the uncertainties given above are reasonable. A fit 
with the center fixed and a solid body function yielded the same systemic velocity. We then 
fixed the center position and central velocity and fit the field in concentric rings that were 
15" in width and stepped every 15". The innermost circle did not yield a solution and that 
point is omitted. For the SW or receding half of the galaxy, we adopted the center position 
and systemic velocity determined from the NE half and again fit in concentric rings. 

The results are shown in Figure [T3] as solid circles. The kinematic axis in each ring 
is traced as the black line on Figure [131 One can see that for each half alone the position 
angle is relatively stable with radius except for the outer two annuli on the SW side. The 
inclination, on the other hand, is quite variable although over much of each half of the field 
it changes systematically. Generally the inclinations are all higher than we expect based on 
the optical major-to-minor axis ratio. The maximum rotation speed is about 15 km s~^ on 
the NE side and a few km s~^ lower on the SW side. 

With a maximum rotation velocity of only ~15 km s^^ (~ 16 km s^^ when corrected 
for an inclination i = 66°) and an average dispersion velocity in the C+Di map of 8.5 km 
s~^ (Figure [151 §4.4p . the galaxy has a significant amount of pressure support from the 
random motions in the gas. This means that simply using the rotation velocity will result 
in an underestimate of the dynamical mass. We have corrected for this by adding in the 
pressure support; this is c alled the asymmetric d rift c orrection. To calculate this, we used 



the procedure described in JBegum &: Chengalun (120041 ) and we follow their discussion here. 



To make this applicable to a gas disk rather than a coUisionless stellar system for which 
the rotational velocities are much larger than the random motions, we make the assumption 
that the pressure support can be approximated by the gas density times the square of the 
velocity dispersion: 

vl = vl - ra^[^(lnSHi) + ^^{\na^) - ^(ln2/^.)]. (1) 

Here, Vc is the corrected circular velocity, vq is the observed rotation velocity from the 
original rotation curve, Shi is the H I surface density, a is the velocity dispersion, and h^ is 
the scale height of the H I disk. We further assume that neither the velocity dispersions nor 
the scale height vary significantly across the disk (i.e. with radius r), so their derivatives are 
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zero, leaving us with 



2 2r 

;„ — rrr \ 
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''dr 



(In Shi)]. 



(2) 



To determine the derivative of the H I surface density, we must determine a functional 
expression to describe it. Towards this end, we fit the profile (derived using the optical 
center) shown in Figure [H] with a Gaussian of the form: 



SHi(r) = Soe-(-'^)'/2^o 



(3) 



where So is the amplitude, /i is the radius at which the Gaussian peaks, and tq is the scale 
length. The profile was fairly well fit with Sq = 11.42 ± 0.54 Mq pc'^, /i = 26 ± 6" (0.55 
kpc), and ro = 57 ± 3" (1.26 kpc). Taking the derivative, eq. [2] now becomes 



Vn + 



(4) 



Using the values found from the gaussian fit and a galaxy-wide average of 8.0 km s^ as a 
conservative lower estimate for the velocity dispersion (corrected for observational broadening 
due to the velocity resolution), we then corrected the rotation curves for both the NE and 
SW halves. The result is shown in Figure [H] (open circles), along with the uncorrected 
rotation curves (solid circles). Applying the correction increases the maximum velocities to 
20 km s~^ for the NE half and 22 km s~^ for the SW half. The corrected rotation curve is 
essentially that of a solid body on both sides of the galaxy. 

The necessity of correcting the rotation to include pressure support is apparent when we 
calculate a dynamical mass for VII Zw 403. Using the highest rotational velocity from the 
uncorrected rotation curves and M = ^ax-^max/C we get an estimated dynamical mass of 
7.9 X 10^ M(7). This is an M}^j/Mc\vn of 0.66, which is hig her than the ranges for BCDs found 
by others ( iTaylor et al.lll994J : Ivan Zee et al.lll998l l200l[ ). The problem becomes worse when 
the total luminous mass (Mgas + M^tars) is considered. The gas mass is found by multiplying 
the H I mass by a factor of 1.34 to correct for He. No correction is made for molecular gas; no 
good estimates exist as CO det ections in BCDs have been few, and thos e primarily in BCDs 
with elliptical inner isophotes ( ISage et al.lll992l ). or in metal-rich dirrs: iTaylor et al.l (Il998[ ) 
report from a literature sample that there are no known detections of CO in galaxies with 
12 + \ok(0/H) < 7.9. VII Zw 403 has irregular inner isophotes and 12 + \og{0/H) = 7.69 
( llzotov &: Thuanlll999[ l. so is unlikely to have detectable CO. 



To find the mass in sta rs, we can estimate a lower limit by using a stellar M/Ly ratio 
from iBell fc de JoneJ feoOlh that is a fui iction of ( ^ - y )n color: log(M/Lv) = -0.37 + 
1.14(i? — V)q. This is appropriate for a ISalpeterl (I19551) stellar initial mass function and 
metallicity Z = 0.008 and is derived using the iBruzual fc CharlotI (|2003[ ) stellar population 
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models. This results in a stellar mass of approximately 3.9 x 10^ Mq, or a stellar mass-to- 
light ratio of 0.89. However, this does not t ake the starburst into account. Assuming that 
the SFR was 30 times higher for 200 Myr (JLynds et al.lll998l ) during the past 13 Gyr, we 
then add that to the constant SFR model (0.013 Mq yr~^) and obtain an upper estimate of 
9 X 10'^ Mq (stellar M/L-q = 1.6). With this range, our luminous mass estimate is between 
1.1 X 10^ and 1.6 x 10^ Mq. Even the lower of these is greater than the highest dynamical 
mass estimated from the uncorrected rotation curve by a factor of 1.4. 

However, when we use the corrected rotation curves to calculate a dynamical mass, 
we get Mdyn = 2.5 x 10® Mq, which is 1.5 times larger than even the higher estimate of 
the luminous mass . Th is is now a M H j/MHyri of 0.21, which is more in line with those in 



Tavlor et al.l fll994J ) and Ivan Zee et al.l (119981, 120011 ) 



4.4. Velocity Dispersions 

Contours from the B+C+Di data showing the velocity dispersions across the galaxy 
overlaid on the H I flux map are shown in Figure [151 In general, the higher dispersion areas 
(~ 10 km s~^) correspond to the area where most of the Ha flux originates, as shown in 
Figure Uni There is nothing notable in either the V or Ha images in the unresolved region 
of highest dispersion (lo cated just NE o f the large H I clump), but it is adjacent to an 
unresolved x-ray source (lOtt et al.ll2005al ). as discussed below. 



Discussion 



5.1. Relation between Star Formation and Gas 



As can be seen in Figure [T71 the highest H I column density region r oughly corresponds t o 
the region of strongest Ha emission. This is seen in other BCDs as well (Ivan Zee et al.lll998l ). 
There is diffuse Ha emission that spreads to the north, but not as far as the H I emission. 
Thuan et al.l (|2004J ) find that there is weak continuum emission (1.4 GHz) associated with 
the Ha emission; with a slight extension to the north at 1.4 GHz as well. Without a spectral 
index, they were unable to verify that the continuum emission is mainly thermal in nature 
(and thus directly due to to the H II region), but it does seem likely. 



Qtt et al.l (l2005al Jbl) presented Chandra x-ray observations of several dwarf starburst 



galaxies, including VII Zw 403. Their data show no diffuse x-ray emission associated with 
VII Zw 403. They found that the diffuse x-ray emission they detected in the other galaxies 
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in their sample was probably associated with the development of galactic winds from coronal 
gas heated by star formation. They suggest that the lack of detected diffuse x-ray emission 
for VII Zw 403 (they give an upper limit of 2.8 x 10^^ erg s~^) could be due to several possible 
factors, including low metallicity or smaller volume. The expanding Ha shells in the two 
galaxies with no detected diffuse x-ray emission (VII Zw 403 and I Zw 18) are 10-17 Myr 
older than those in the galaxies for which they do detect diffuse emission. Perhaps a coronal 
wind never developed, or perhaps it has already had time to cool. 



Ott et al.l ( l2005a[ ) did detect one unresolved x-ray source whose power-law spectral index 



is consistent with that of an x-ray binary, but did not detect the filamen ts or lobes seen by 



Papaderos et al.l (119941 ) from ROSAT data. Similar x-ray sources in the lOtt et al.l (l2005af ) 
sample of 8 starbursting galaxies are located near "bright H II regions, rims of superbubbles 
or young stellar clusters." The location of the unresolved x-ray source relative to the Ha 
emission can be seen in Figure [T71 The x-ray source is located near, but not coincident with, 
an Ha knot. As seen in Figures [16] and [TTl the x-ray source is also on the north edge of the 
high dispersion knot, and is located in/on the northeast side of the depression/cavity in the 
H I, in a smaller region of even lower density. Recall that the broadband surface photometry 
(Figure E]) was consistent with constant star formation in the outer regions over the past 
10-20 Gyr, with more recent star formation (past 200 Myr) in the central regions. If the 
x-ray source is, indeed, an x-ray binary, then perhaps we are seeing an evolved stellar system 
located near a region of recent star formation, and also near a region of disturbed neutral 
gas. The logical explanation is that a star formation episode and/or resulting supernova has 
injected energy into the H I, increasing the dispersion; and depleting or perhaps clearing out 
some of the neutral gas. 

To test this, we examined the kinematics of the gas in the B+C+Dq cube in a small 
region around the x-ray source. To enhance the signal-to-noise in the spectra, we "collapsed" 
a six-pixel square region centered on the x-ray source to a single pixel, and then examined 
the resulting spectrum. There was a fairly clear double-peak visible, so we fit a pair of 
gaussians to them (using XGAUS in AlPs). The two peaks are separated by 13.5 km s^^. 
This is consistent with expansion (or contraction) of ~ 7 km s~^. The apparent radius of 
the depression around the x-ray source in the B+C+Dq data is ~ 105 pc, which is roughly 
the same size as the beam. This means that the region may not be resolved, so this size is 
an upper limit. If we assume a constant expansion of 7 km s~^, we get approximately 15 
Myr to clear out a depression of that size. 

We then t ried to estimate t he amount of energy required to form such a depression using 



the method in IChevalierl (119741 ). Here we assume a spherically symmetric expansion. If the 



expansion has already blown out of the disk, this may not be the case (and would increase 
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the estimated age of the feature; see below). For a velocity of 7 km s~^ and a radius of 105 
pc, we find the total energy required to be Eg = LQuq^"^ x 10^^ ergs. To get an estimate for 
the density, uq, we assumed a pre-existing column density of 1.5 x 10^^ cm~^ based on the 
average density in the region around the small depression, and a scale height for the galaxy 
of 200 pc. The scale height was taken to be approximately equal to the radius of the larger 
H I depression/cavity (see below). This gives us no = 1.0 cm~^, and a resulting total energy 
of 1.6 X 10^^ ergs. This is similar to the 10^^ ergs from a typical Type II supernova, s o it is 



feasible in this sense. We then attempted to apply the models of iMcCray fc Kafatod (119871 ) 
for the formation of supershells from stellar winds and supernovae to estimate the number of 
stars with masses greater than 7Mq required to produce a hole of a given size. The resulting 
number was approximately 200, of which about 50 should be O stars. However, the resulting 
timescale for shell formation calculated from the models was short, only 2 Myr. In this short 
time, the expansion due to the star formation is still primarily from stellar winds; and it is 
unlikely that any stars would yet have exploded as supernovae, and we should be able to see 
an OB association of this size — which we don't. So about all we can conclude is that the 
x-ray source appears to be coincident with an expanding depression in the H I; and adjacent 
to a region of disturbed gas to the south. 

We also examined the kinematics of the gas around the edges of the larger depres- 
sion/cavity in the H I. However, no expansion was detected. If we now use the average 
velocity dispersion in the gas (9 km s~^) as an upper limit to any expansion, with a radius 
of 205 pc we obtain a lower limit of 22 Myr for a hole this size to form. Because we don't 
detect expansion, the hole could be stalled because it has broken out and is freely expanding 
perpendicular to the disk, and so could be older than that. We used the same methods to 
estimate the required energy and number of stars required to form such a depression as for 
the area around the x-ray source. This time, we used a column density of 2.5 x 10^^ cm~^ 
based on the existing gas density around the cavity; giving us an no of 1.6 cm~^. The re- 
quired energy is the r i 3 x 1 0^^ ergs (equivalent to approximately 30 SNe). However, from the 



McCray fc KafatosI (119871 ) equations, the estimated number of stars needed to create such a 
large depression via winds (and possibly supernovae, but see the timescale discussion, below) 
with masses greater than 7Mq is then 2800, which is extremely high. If there were a young 
cluster that large, it would be easily visible. This suggests an older age for the formation of 
the large H I cavity; or perhaps the gas was slowly depleted/blown-out over long time-scales 
by a series of star-forming events rather than by one large cluster. 



A more serious problem with applying the lMcCray fc KafatosI (119871 ) model is that doing 



so results in a calculated timescale of only 3 Myr; which is not only short in stellar evolution 



terms, but also inconsist ent with the 



ower limit of 22 Myr we calculated, above. A similar 



disconnect was found by lLynds et al.l (119981 ) for formation of the large H II regions in VII Zw 
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403. They used the Wide Field Planetary Camera 2 (WFPC2) on the Hubble Space Telescope 
(HST) to produce color-magnitude diagrams (CMDs) of several of the bright H II regions 
in VII Zw 403. The brightest H II region (visible as the darkest region of the grayscale 
in Figure [T7] and marked in Figure [TK]) was consistent with an OB association containing 
approximately 25 bright blue (O type) stars; the age of the cluster based on the CMD was 
4-5 Myr. From their ground-based long-slit emission line spectra of the H II regions, they 
were able to detect expanding Ha shells; including a partial shell around the brightest H II 
region with a diamet er of 79 pc and an expan sion velocity of ~ 65 km s^^. They then 
attempted to use the iMcCray &: Kafatoj ( 119871 ) models as we did, but they also found a 
timescale that was too short: ~ 1 Myr, which did not match the 4-5 Myr age estimate from 
the CMD. They concluded that p erhaps this indicated that the starburst that formed the 
association was not instantaneous ( IShuU fc Sakenlll995[ ). 



Lynds et al.l (119981 ) also produced CMDs of some of the central areas of VII Zw 403. The 



location of their "center" and "north" fields are marked in Figure [181 which shows the HST 
image with contours from the B+C+D_i integrated flux map. Neither field corresponds 
exactly with the H I cavity/depression, but the "north" field is located on the east side 
of it, and just south of the x-ray source. The CMD for the "north" field is less rich in 
massive stars than the one for the "center," which revealed approximately 20 red supergiants 
(> 10 Myr) and some blue supergiants (~ 5 Myr), but both regions cont ain a mixed-age 
population of stars with ages of 5 and 10 Myr. In general, iLynds et al.l (Il998[ ) find a relatively 
smooth elliptical distribution of the evolved stars with a stellar surface density that falls off 
exponentially. The younger, bluer massive stars tend to exhibit clustering into associations 
and groups, and are primarily located in the central regions of the galaxy. The youngest 
are associated with the bright H II regions, which we find are co-located with H I peaks 
(Figure [T9|) . They posit star formation 1-2 Gyr ago followed by a strong starburst episode 
600-800 Myr ago, with less vigorous star-formation since then. Interestingly, the kinematic 
center of VII Zw 403 as identified by the rotation curve analysis lies close to the center of 
the large H I cavity. If the cavity is the result of star-formation, this would be consistent 
with centrally-concentrated star-formation activity in the past. There is evidence that stellar 
feedback from star forr nation events can an d does create sometimes large (~ 1 kpc) H I holes 
and shells, as found by IWeisz et all J2009I ) in IC 2574. 



Another possibility for forming the hole is consumption of the gas by star formation. 
The timescale to consume the gas that filled the hole at the current star formation rate, 
given our estimate of the pre-existing level of the gas in that region, is of order 700 Myr. 
This is relatively short if star formation has been on-going in this region for that extended 
period of time. Alternatively, the original gas in this region could have been co nsumed quite 
quickly by the starburst that took place 600-800 Myr ago (JLynds et al.lll998[ ). but we then 
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require that the hole not be replenished very rapidly, which seems unreasonable given that 
the hole should have filled in long ago just from the dispersion in the gas. 



5.2. Origin of the Starburst 

One possible exp lanation for t he en hanced star formation episode that characterizes 
BCDs is proposed by ISalzer et al.l ( Il999l ). who suggest that it is a natural consequence of 
the central mass density concentration found in BCDs. In an optical study of a sample of 18 
BCDs and 11 dirrs, they find that when the luminosity contribution from the starburst is 
ignored (by fitting the underlying host galaxy's surface brightness profile as determined from 
an exponential fit to only the outer portions), the underlying hosts of BCDs have shorter scale 
lengths and higher central surface brightnesses than dIrrs with similar absolute magnitudes. 
They cite the shorter scale lengths as evidence of stellar compactness: that BCDs are low- 
surface brightness galaxies with the highest central stellar mass concentrations in the dwarf 
spectrum. 

There is s ome e vidence that this may be the case for the H I gas in BCDs as well: 
van Zee et al.l (119981 ) plotted H I surface density profiles for a sample of BCDs and dIrrs, 
and found that the gas in BCDs is more centrally concentrated than in similar dIrrs. It is 
commonly thought t hat the gas surfa c e density mu s t exce ed a critical threshold in order to 
support a starburslo ( lKennicuttlll989l ): ISalzer et al.l (119991 ) also suggest that the high central 
mass density conce ntration allows these galaxi es to retain their gas, allowing episodic bursts, 
and calculations by lMac Low &: Ferraral (1l990l ) indicate that mass loss should be minimal for 
galaxies with gas masses gr eater than 10^ Mq. A follow-up study of the H I in a sample of 
BCDs by lSalzer et al.l (120021 ) found that the Mhi/Lb for these galaxies is about twice as high 
as that for low surface brightness dwarfs when only the luminosity of the underlying low- 
surface brightness component of the BCD is taken into account. This provides a significant 
gas reservoir to support the starburst episodes, and could enhance the ability of the galaxy 
to retain gas during/after a burst. Additionally, if the gas remains concentrated even after 
a burst, then the timescale between bursts could be short. This would make it difficult to 
find a BCD in the "off" state. 

If this hypothesis is correct then we would expect VII Zw 403 to exhibit signs of cen- 
tral compactness, both optically and in the H I. We find that compared to the sample in 



^Although azimuthal averages of the gas density in dwarfs ra rely exceed this cr i tical value, there is 



evidence that local values in the region of ac t ive star formation do (Hunter et al.ll200ll : Ivan der Hulst et al 
199i Ivan Zee et aDll997t JMeurer et aLlll998t Ide Blok fc Walteijbooei ). 
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Salzer et al.l ( 119991 ). VII Zw 403 has a slightly lower fiso at a given M^ than their sample 
of BCDs, although slightly higher than for their LSBDs. As mentioned in fj3l VII Zw 403's 
scale length (relative t o its My) lies at the peak o f the distribution of the sample of dims 
and BCDs plotted in iHunter &: Elmegreenl ( 120061 ). Compared to the ISalzer et al.l (Il999[ ) 
sample, however, it is large not only for their sample of BCDs, but also compared to their 
LSBD sample. So the underlying stellar component does not appear to be greatly centrally 
concentrated. However, there are some indications for central concentration in the H I on 
larger scales (i.e. at a linear resolution similar to that for which structure is seen in the H I 
in dirrs; see Figure |9]). 

However, since the vigorous star formation episodes that characterize BCDs are rela- 
tively short-lived compared to the age of the host galaxy, p e rhaps these systems don't start 
out with a high central mass concentration. IVan Zee et al.l ( 120011 ) finds that BCDs are low 
angular momentum systems, and suggests that they have collapsed to more centrally con- 
centrated systems with steeper rotation curves. This affects the critical threshold density 
for star formation, and makes BCDs more susceptible to star bursts than higher angular 
momentum systems. This "collapse" requires a redistribution of mass in the galaxy. One 
proposed mechanism for mass redistribution and angular momentum transfer in low mass 
BCDs (without nearby companions) is a barlike torque arising from the irregular symmetry 
and i rregular mass distribution in these small systems, or from an elongated dark matter 
halo ( JHunter fc Elmegreenl 120041 ). Such mass flow should show up in the rotation pattern of 
the gas disk as a central twisting of the isovelocity contours. VII Zw 403 has a change of 
15-20° between the NE and SW parts of the galaxy, so we can't rule out the possibility of 
some kind of mass flow in the gas disk. 

Another possible mechanism for flattening the gravitatio nal potential and expansion of 
the stellar component is suggested by lPapaderos et al.l (Il996bl ). who propose that the heating 
of the ISM by the energy from the starburst could be the mechanism for mass and angular 
momentum redistribution. Although VII Zw 403's scale length is not particularly short, 
which could indicate that it has undergone some expansion, the lack of diffuse x-ray as well 
as the lack of extensive large-scale structure in the H I (the H I depression discussed above 
occupies only a small, central part of the galaxy) would seem to indicate that not much 
energy has been injected into the galaxy as a whole by star formation. It seems unlikely, 
then, that enough energy would have been produced to initiate any signiflcant redistribution 
of the mass in the galaxy, especially the stellar component. In conclusion, it isn't clear that 
this theory of large-scale mass re- arrangement can explain what we see in VII Zw 403. 

Interactions, either with other galaxies or with small gas companions, h ave also been 



proposed as the triggering mechanism for the starbursts that occur in BCDs (see lBrosch et al. 



-30- 



2004 for an overview of the literature), yet surveys to detect nearby possible perturbers 



around BCDs have had mixed resu 
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Hunter fc Elmegreenl |2004J ). VII Zw 403 would seem to be in the set 



of galaxies that haven't been tidally triggered, as we have found no nearby H I clouds or 
companion galaxies. However, if it had accreted a small perturber in the distant past, this 
might be difficult to detect. 

There is observational evidence that some dwarf systems have under gone accretion or 
merge r events in the past. A nearby example is thought to be IC 10, which IWilcots fc Miller 
( 119981 ) argue still has primordial gas settling into the galaxy via accretion. Classic signatures 
of accretion include the appearance of two distinctly different kinematic components, visible 
in the H I or stellar velocity fields. Although the gas isovelocity contours for VII Zw 403 are 
somewhat perturbed, this isn't unusual for BCDs ( §4.3p . and there don't appear to be two 
completely separate components despite the change in position angle of the H I kinematic 
axis. Another ac c retion signature is a misalignment between the optical and H I axes. 
Hunter fc Wilcotd (120021 ) plotted the difference b etween optical and H I kinematic axes for 
47 Im and Sm galaxies taken from ISwaterd (119991 ) and found that about half had offsets less 
than 10°, with the maximum being 40°. For VII Zw 403, the optical axis is at —11°. In the 
NE half of the galaxy, the H I kinematical axis is at +45° to +50°, and is around +20° to 
+30° in the SW part, so there is a somewhat significant misalignment between the gas and 
optical components of between 30° and 60°. 



For further comparison, we look to the gas-rich Im galaxy DDO 26, for which iHunter fc Wilcots 



(120021 ) conclude that the most likely explanation for the lack of relaxation they detect in its 
gas kinematics is the on-going merger of two small gas-rich objects. They find an optical/H I 
axis misalignment of 60°, along with a change in the H I position angle with radius. They 
also find that the velocity field of DDO 26 is fairly regular in the inner regions and becomes 
less so in the outer disk. These are features that are similar to but more extreme than those 
we see in VII Zw 403. Differences include the presence of a faint H I arm in DDO 26 and the 
detection of two separate peaks in the H I intensity which were both fairly broad (~ 40 km 
s~^) in velocity. Additionally, DDO 26 is currently quiescent, with a SFR of 0.0004 Mq yr~^ 
kpc~^ compared to VII Zw 403's rate of 0.02 Mq yr~^ kpc~^. So from these comparisons, it 
seems possible that VII Zw 403 underwent an accretion or merger event in the past. 

Another possibility to explain the skew in the H I kinematic axis is that expansion of 
the H II regions has caused a subsequent expansion in the surrounding H I, lifting it out of 
the plane of the galaxy. If the near side of the galaxy is to the west (right), and we envision 
the expanding H I as distributed in a torus aligned with the plane of the galaxy, then gas 
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to the left of the kinematic center would be blue-shifted relative to the bulk of the galaxy. 
Likewise, the gas to the right of center would be somewhat behind the bulk of the galaxy and 
red-shifted. The expanding gas would then skew the observed velocity field in the sense that 
we see: the isovelocity contours in the north part will be skewed to the left, and those in the 
south skewed to the right. If the expansion is on the order of the 7 km s~^ we measured for 
the lowest density region around the x-ray source, or the 9 km s~^ limit set by the average 
dispersion, then in this slowly rotating galaxy, the expected skew would be large, which is 
what we see. In this case, the hole in the H I would be due to a combination of consumption 
and slow blow-out. 



6. Summary 

We have used multi-wavelength data to examine the relation between the stars and 
gas in VII Zw 403 in an effort to understand the origin of the starbursting activity that 
characterizes BCD galaxies. As for many BCDs, in the optical V^-band we find a mostly 
symmetrical, elliptical distribution of stars; this is overlaid with a color gradient that is 
bluer in the central regions and redder in the outer. This indicates that the central regions 
have undergone star formation in the past 200 Myr, whereas the outer region colors indicate 
constant star formation over the past 10-20 Gyr. The galaxy also has a few large H II 
regions concentrated toward the center of the galaxy (within the inner 1.3 kpc) where star 
formation has occurred in the past 10 Myr. Most of the activity is in a large complex just 
south of the central part of the galaxy. Overall, the SFR for VII Zw 403 is at the high end 
of the normal range for dirrs, which is, again, typical for BCDs. However, the galaxy can 
continue to form stars at this rate for at least 5 Gyr more; 10 Gyr if we include gas recycling 



( JBrinchmann et al.ll2004J ). 



Our new VLA data in combination with earlier data resulted in a lower total H I flux 
detection than single-dish estimates; we believe this is due to inclusion of Milky Way emission 
in the single-dish spectra and in the low-resolution VLA observations. We have attempted 
to remove this Galactic contamination from our low-resolution (D configuration) data for 
our analysis here. Our total H I mass estimate then is 5 x 10^ Mq, which is lower than 



the median for the Salzer BCD sample (jSalzer et al.l 120021 ). but within the normal range of 



values. However, the VIhi/^b indicates that VII Zw 403 is somewhat gas-rich for a BCD of 
similar luminosity. Co nversely, it rnay be "luminosity-poor" relative to other similarly-sized 



BCDs: earlier work by lLynds et al.l (Il998[ ) concluded that the SFR for VII Zw 403 was up to 
30 times greater 600-800 Myr ago. Perhaps we are seeing it post-burst. The idea that VII 
Zw 403 is presently past the peak of a starburst event could also explain the non-detection 
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of diffuse x-ray emission by lOtt et al.l ( 2005a| ]bl). indicating a lack of coronal wind or that it 
has had time to cool. 

The H I velocity field, although complex with a change in position angle of the kinematic 
axis between the NE and SW components of about 15-20°, is not that unusual for a BCD. 
We see no other strong indications of a large-scale gravitational interaction in the velocity 
field, such as tidal tails or strewn-about H I. The skew in the velocity field can, however, be 
explained if the H I in the galaxy is undergoing a general symmetric expansion away from 
the plane on both sides surrounding the star-forming region as a result of pressure from the 
Ha region. The rotation curve is essentially solid-body for both halves of the galaxy, and 
with correction for as ymmetric drift, produces a dynamical mass e stimate that is similar to 
those of other BCDs JTavlor et al.lll994l : Ivan Zee et aPllQQsl . boOlh . 



If the galaxy is post-starburst, then one might expect to see some effect in the inter- 
stellar medium from the energy injected by stellar winds and supernovae. Additionally, star 
formation should use up some of the neutral gas in the galaxy. In accordance with this, 
although the H I is more centrally concentrated in VII Zw 403 than in d i rrs, th e central H I 
surface density is lower than is typical for the BCDs in the Ivan Zee et al.l (119981 ) sample, and 
the H I is more extended — out to 4.3 R25 (2.3 R^). At large scales (650 pc), the H I distribu- 
tion appears relatively smooth as for most BCDs. This is in contrast to dirr galaxies, which 
at the same scale often show structure in the gas component in the form of holes, shells, and 
ridges. For VII Zw 403 however, we begin to see some structure in the H I distribution only 
at small (200 pc) scales. 

There is clearly a connection between recent star formation and the H I distribution and 
kinematics in the galaxy. The youngest stars (4-5 Myr old) are located near H II regions 
and local peaks in the H I, with the regions of highest H I column density corresponding 
to the regions of strongest Ha flux (as seen in other BCDs). Our H I maps also show 
that the regions of higher dispersion are located near, although not completely coincident 
with, with Ha regions; perhaps in response to the energy injected by the star formation 
activity. Interestingly, the region of highest dispersion in our map is directly adjacent to the 
unresolved x-ray source detected by lOtt et al.l (l2005al ) that could be an x-ray binary. The 
source is located near an H II region, and is in a small lower-density region at the edge of a 
larger depression or cavity in the H I. The center of the large H I cavity is coincident with 
the kinematic center of the H I velocity field. 

We attempted to investigate whether the small low-density region could be the result 
of star formation/SNe blowing a hole in the gas. Although both the amount of energy 
required and the number of massive stars needed were reasonable, the timescale from the 
shell formation model was too short for many SNe to have already occurred. We performed 



-33- 



the same analysis for the (stalled) large cavity in the H I, but found that in addition to 
requiring the presence of a large star cluster (not seen), the ti mescale froni themodel was 



again too short. Similar timescale discrepancies were found by iLynds et al.l (Il998l ). leading 



them to invoke non-instantaneous starburst/star-formation activity in a series of events. 

So we are brought back to the question of what triggered the vigorous star formation 
episode that apparently peaked in the past but is still ongoing today. We have investigated 
three possible triggers: tidal forces from an external perturber, inherent central mass density 
concentration, and accretion of a smaller object. We found no nearby large or small per- 
turbers, so it has apparently not been tidally triggered. VII Zw 403's mass doesn't appear 
to be strongly centrally concentrated either, as evidenced by its optical scale length, lower 
central surface brightness, and lower H I surface density than comparable BCDs. This leaves 
us with an accretion event. There is evidence for accretion in some dwarf systems, including 
some Local Group dwarfs. Compared to another suspected accretion/merger system (DDO 
26), VII Zw 403 has some similarities but less extreme signatures. Of the three possible 
triggers we have investigated, this seems the most likely; although the evidence is far from 
conclusive. 



This publication makes use of data products from the Two Micron All Sky Survey, 
which is a joint project of the University of Massachusetts and the Infrared Processing and 
Analysis Center/ California Institute of Technology, funded by the National Aeronautics and 
Space Administration and the National Science Foundation. This work is also based in part 
on archival data obtained with the Spitzer Space Telescope, which is operated by the Jet 
Propulsion Laboratory, California Institute of Technology under a contract with NASA. 

We appreciate the use of Ohio State University's Fabry-Perot and the Ohio State In- 
frared Imager-Spectrograph and Mark Wagner for making them work for us. The authors 
would like to thank the anonymous referee, who provided valuable comments and criticism 
that improved the paper. CES would like to thank Glenn Morrison and Frazer Owen for 
help with data calibration questions. 

Support for CES came from grants AST-0407051 and AST-0707468 from the National 
Science Foundation. Support to DAH and TEN for this research came from the Lowell 
Research Fund and grants AST-9802193 and AST-0204922 from the National Science Foun- 
dation. Support to DAH also came from NASA through grant NNX07AJ36G. Funding to 
BGE was provided by NSF grant AST-0707426. Support for TA came from the National 
Science Foundation Research Experience for Undergraduates program (AST-0552798) in 
partnership with the Department of Defense ASSURE programs (Awards to Stimulate and 
Support Undergraduate Research Experiences). 



-34- 



Facilities: VLA GALEX Lowell Observatory 
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Table 1. VLA Observations 



Parameter 




D configuration 


C configuration 


B configuration 


Observation Date 




1997 Nov. 10 


1992 April 


11 


2006 Sept. 10 


Time on Source (min) 




253 


220 




546 


Central Velocity (km s^^) 




-60 


-60 




-60^ 


Bandwidth (MHz) 




1.56 


1.56 




1.56 


No. of Channels 




128 


128 




128 


Velocity Resolution (km s^ 


') 


2.6 


2.6 




2.6 


Beam Size (arcsec) 




72.91 X 53.57 








Single Channel r.m.s. (mJy 


beam~^) 


0.93 








Single Channel r.m.s. (K) 




0.28 









^Doppler tracking not available; central fixed frequency chosen to correspond to —60 km s ^. 
See text for details. 



^Natural weighting 
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Table 2. Combined Configuration Parameters 



Parameter C+Di B+C+Di B+C+Dq B+C+D_i 

Robustness Factor +1 +1 —1 

Beam Size (arcsec) 34.10 x 27.60 12.85 x 11.71 9.23 x 7.57 4.95 x 4.63 

Beam Size (pc) 730 x 590 280 x 250 200 x 210 110 x 990 

Single Channel r.m.s. (mJy beam^i) 0.83 0.80 0.86 0.90 



-37- 



Table 3. Summary of Properties. 



Parameter 



Value 



R.A. (B1950; V) 11'^ 24"^ 33*^4 ±4" 

Decl. (B1950; V) 79° 16' 10"± 3" 

R.A. (B1950; HI) ll'^ 24^" 35^1± 10^3 

Decl. (B1950; HI) 79° 16' 18"± 27" 

b/a (V) 0.49 

P.A. (V) -11° 

D (Mpc) 4.4 

Mhi (Mq) 5.23 X 10^ 

Rui'^ (arcsec, kpc) 157, 3.35 

E{B-V)f'° 0.023 

R25 (arcsec, kpc) 36.5, 0.78 

Rh (arcsec, kpc) 68.95, 1.47 

fj,Q (magnitudes arcsec^^) . 23.12±0.09 

R^ (kpc) 0.50±0.01 

Mvo (magnitudes) -14.26±0.009 

(U -B)o -0.50±0.01 

(B -V)o 0.27±0.01 

(\/-J)o 1.30±0.01 

(FUV-NUV)o 0.17±0.01 

(NUV-V)o' 0.70±0.01 

log LHa,o (ergs s^^) 39.34±0.001 

log Lnuv,o (ergs s^^ Hz~i) 26.07±0.001 

SFR^/ ' (Mo yr-i) 0.015 

SFRi^c// (Moyr-i) 0.013 

SFRcHa'^ (Mo yr-i kpc-2) -1.72 



^Measured to A^hi 



1 X 10^9 cm-2. 



^E{B — V)f is foreground reddening due to the 
Milky Way (Burstein & Heiles 1984). For the stars 
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in VII Zw 403, we assume an additional internal red- 
dening of 0.05 magnitude; for the H II regions we 
assume an additional internal reddening of 0.1 mag- 
nitude. 

'^NUV is an AB magnitude. 

'^Star formation rates derived from L Hn or t he ul- 
traviolet using the formulae of iHunter et al.l (l2010f ) 



that integrates from 0.1 Mq to 100 Mq with a 
Salpeterl (jl955[ ) stellar initial mass function. 



'^Star formation rate normalized to the area within 
one scale-length R^. 
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Fig. 1. — Images of VII Zw 403 in various passbands. The images are shown as the logarithm 
in order to allow the inner and outer structure to be viewed simultaneously. The images are 
also geometrically transformed to match the scale and orientation of the V image, and the 
same field of view is shown in each plot. North is up and East to the left. 
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Fig. 2. — Single-dish spectrum from lThuan fc Martini ( 1l98ll ). 
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Fig. 3. — Summed flux in each channel from the uniformly- weighted D configuration data. 
The strong spike and absorption features near km s~^, as well as some of the emission from 
approximately —55 to —80 km s^^ are due to Galactic H I. 
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Fig. 4. — Summed flux in each channel from the naturall y-weighted D configurat ion data 
(solid line) is plotted along with the single-dish detection by lThuan fc Martini ( 1l98ll ) (dotted 
line). The emission located around —70 km s~^ is foreground H I from the Milky Way. 
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Fig. 5. — Channel maps from part of the D configuration data cube. Contaminating emission 
from the Milky Way is visible in the velocity range from —23.9 to —80.6 km s~^. The strong 
concentrated emission in the center of the field is from VII Zw 403. 
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Fig. 6.— Azimuthally-averaged V, Ha, NUV, 3.6 /im, {U - B)o, {B-V)o, {FUV-NUV)o, 
and {V — J)o surface photometry and H I surface density of VII Zw 403, using the V-band 
morphological center of the galaxy. The UBVJ and NUV photometry are corrected for 
reddening as described in the text. The exponential fits to the V , 3.6 /im, and NUV-band 
surface brightness profiles are shown as solid lines in the top panel. These profiles are fit with 
two exponentials; the outer one shallower than the inner. The Ha luminosity is proportional 
to the star formation rate. The scales for T^ua and for Hhi have been set to cover the same 
logarithmic interval as the broad-band surface photometry covers in magnitudes and are 
shown on the right y-axis. 
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Fig. 7. — Star formation rate per unit area (S FRn; from iHunter fc Elmegreenl ( l2004j ) for 
unbarred dim galaxies. After Figure 11 from iHunter fc Elmegreenl (120061 ). SFR/) is the 
integrated Ha-based star formation rate divided by the area in one l^-band disk scale length 
of the galaxy. VII Zw 403 has a star formation rate that is significantly higher than the 
median in this sample. 
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Fig. 8. — H I surface density in VII Zw 403 from the C+Di data. The line with the solid 
circles labelled "optical" used the center, inclination (66°), and position angle (—11°) from 
the V-band data. The line with open circles labelled "H I" used the center, inclination 
(77°), and average position angle (47°) from the H I rotation curve for the NE half of the 
galaxy. H I surface density profiles for galaxies with simila r beam/-Ri/3 or beam /-R25 ratios 
are also sho wn. Data for these profi les are from DDO 43 (jSimpson et al.ll2005bl) : DDO 75 
fSextans A) Jwilcots fc Huntei] l2002h: DDO 154 JCarignan fc Beaulieulll989h : and DDO 155 



(jCarignan et al. 



19901). 
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Fig. 9. — The V image overlaid with contours from the C+Di configuration integrated fiux 
map. The beam size is indicated in the lower left corner. The levels are (2, 10, 20, 30, 40, 
50, 100, 150, and 200) xlO^^ cm~^. The lowest contour is at l.Scr. (A single-channel 2a 
detection integrated over 3 channels is equivalent to a column density of 1.5 xlO^^ cm^^, 
which we take as indicative of the noise level in the integrated map.) The field of view is 
approximately 5'x 5'. 
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Fig. 10. — Integrated flux density map from the B+C+Dq configuration data. Tlie beam 
size is indicated in the lower left corner. The contour levels are (5, 10, 15, 23, 30, 40, 50) 
xlO^° cm~^. The lowest contour is at 2.4o". (A single-channel 2a detection integrated over 3 
channels is equivalent to a column density of 2.12 xlO^*^ cm^^, which we take as indicative 
of the noise level in the integrated map.) The field of view is 1.9' x 2.3', which corresponds 
to 2.4 X 2.9 kpc. 
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Fig. 11. — Channel maps of the H I emission in the C+Di configuration data cube. The 
velocity of each channel is printed along the top; the beam size is indicated in the first (top 
left) channel. The contours are at (3, 10, 20, 40, 60, 80, and 100) x the single-channel r.m.s. 
(0.83 mJy beam~^). 
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Fig. 12. — H I velocity field from the C+Di configuration data. The beam size is indicated 
in the lower left corner. Contours are in km s^^. 
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Fig. 13. — Isovelocity contours from the C+Di configuration data moment one map su- 
perposed on a V-band image of VII Zw 403. The contours, in white, are labeled with the 
velocity in km s~^. The black line running more or less perpendicular to the isovelocity 
contours shows the position angle of the rotation axis in 15" steps. The plus sign marks the 
kinematic center of the galaxy, and the open circle marks the center in the V-band. Note 
that the uncertainties on the kinematic center {^ 27") are approximately the same size as teh 
FWHM of the beam, which is indicated by the ellipse in the lower left corner (34.1"x27.6"). 
The velocity resolution is 2.6 km s"^. 
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Fig. 14. — Best-fit rotation curve for VII Zw 403. The velocity field was fit in two separate 
parts: the NE or approaching half and the SW or receding half. Only the central position 
and systemic velocity were fixed in these fits, and the rotation speed, position angle, and 
inclination were determined in 15" wide annuli stepped in 15" intervals from the center. The 
solid symbols are the original rotation curve; the open symbols show the curve corrected for 
pressure support (asymmetric drift correction). 
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Fig. 15. — Velocity dispersion contours superposed on the integrated flux from B+C+Di 
data. The beam size is indicated in the lower left corner. Contours are 4, 6, 8, 10, 12, and 

14 km s~^. 
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Fig. 16. — The Ha image with contours from the B+C+Di velocity dispersion map. The 
beam size is indicated in the lower left corner. Contour levels are 4, 6, 8, 10, 12, and 14 km 
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Fig. 17. — The Ha image overlaid with contours from the B+C+Dq configuration data. 
The beam size is indicated in the lo wer left corner. Th e plus sign marks the location of the 
unresolved x-ray source detected by lOtt et al.l ( l2005al ) . The contour levels are the same as 
for Figure [TOl 
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Fig. 18. — The B+ C+D_i H I integrated flux contours on the HST WFPC2 image from 



Lynds et al.l (119981 ). The large circles show the areas for which they determined their north 
and central region CMDs. The small c i rcle la belled "1" shows the location of the brightest 
H II region examined by iLynds et al.l (119981). The plu s sign indicates the location of the 
unresolved x-ray source detected by lOtt et al.l (l2005al ): the x marks the location of the 
kinematic center as determined from the C+Di data. The beam size, shown in the lower left 
box, is 110 X 100 pc. Contours are at (5, 10,20,30,40,50) x 10^° cm"^. The lowest contour 
is at 0.7a, the second at 1.5cr. (A single- channel 2a detection integrated over 3 channels is 
equivalent to a column density of 6.76 xlO^" cm~^, which we take as indicative of the noise 
level in the integrated map.) 
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Fig. 19. — The B+ C+D-i H I in t egrat ed flux contours on our Ha image. The circles show 
the areas for which iLynds et al.l (Il998[ ) determined their north and central regi on CMDs 
The pl us sign indicates the location of the unresolved x-ray source detected by lOtt et al. 
( l2005al ): the x marks the location of the kinematic center as determined from the C+Di 
data. The beam size, shown in the lower left box, is 110 x 100 pc. Contours are the same as 
for Figure [181 



